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ABSTRACT 


This is the continuation of the development of a computer 
based radome design tool. Measured electric field patterns 
for a AGM-88 High Speed Antiradiation Missile radome were used 
to validate the computer model. This computer model is based 
on a method of moments solution of the E-field integral 
equation for bodies of revolution. 
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I. IHTRODUCTION 


A radome is a dome-like structure designed to provide 
environmental protection for an antenna. Missile radomes in 
particular must be designed to meet the stress of high speed 
flight maneuvers and aerodyncunic heating. At the same time 
they must be able to provide suitedDle electrical 
characteristics for the antenna to receive and transmit 
electromagnetic waves. Nowhere is this more important than in 
the design of anti-radiation missile radomes. These missiles 
must counter a wide frequency spectrum of threats, thus 
requiring radomes with broadband electrical characteristics. 
However, due to the severe structural requirements, the 
electrical characteristics generally become a secondary 
priority in the design. 

Most radomes for missile applications are based wholly or 
in part on the ogive shape. The electrical design of a curved 
radome is complicated because of the many possible 
transmission paths through the radome due to reflections. Ray 
tracing, as illustrated in Figure 1.1, is commonly used to 
analyze radome performance, but this method becomes tedious 
and inaccurate when higher-order reflections are present. Ray 
tracing models also assume a radome curvature large enough so 
that the surface can be considered locally flat. Furthermore, 
ray tracing techniques assume the radome is in the far field 
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of the antenna which Is not necessarily the case in practice. 
As a result, ray tracing techniques are not the most accurate 
models to use in analyzing the performance of an antenna with 
a radome. 


ray tracing techniques 
assiune rays are 
parallel and that the 
scattering surface 
is locally flat. 

defocussing 



vary rapidly in 
direction. 


Figure 1.1 Radome ray tracing approximations. [After Ref.l] 
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The method of moments (fM) is a solution to an integral 
equation that is rigorous in principle. It includes multiple 
reflections, surface waves, and can handle radomes in the 
near-field of the antenna. A MM computer code was developed 
by Francis [Ref. 1] and improved by Klopp [Ref. 2 ]. The code 
was verified for simple radome shapes with single-layer walls, 
but has not been used to model multilayer walls. It is the 
objective of this thesis to con 5 >are measured and computed 
results for a multi-layer radome designed for use by a modern 
anti-radiation missile (ARM) . Verification of the code with 
measured data will allow it to be used with confidence in the 
design process, with a potential savings of time and money. 

Chapter II discusses the history and development of the 
ARM and includes a brief explanation of its seeker and the 
required radome characteristics. 

Chapter III explains the basic elements of radome design. 
Aerodyncimic considerations will be discussed and the various 
tradeoffs that must be considered in the design are addressed. 
The geometry of an ogive radome will also be illustrated and 
explained. 

In Chapter IV, a description of the antenna measurement 
facility, the test procedure, and measured VSWR using the HP 
8510 network analyzer will be discussed. 

Concurrent with the hardware testing, computer modeling 
was performed using the simulation developed by Francis [Ref. 
1] and Klopp [Ref. 2]. The theory behind this simulation is 
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derived in Chapter V. This simulation provides pattern and 
gain estimates as a function of radome parameters and 
geometry. It is based on a method of moments solution of the 
electric field integral equation for bodies of revolution (BOR) 
developed by Mautz and Harrington [Ref. 3] . An important 
capability of the code is that the radome can be located in 
the near field of the antenna. 

Chapter VI presents the measured data from both the 
hardware measurements and the computer simulation. Finally 
Chapter VII will present some conclusions and guidelines for 
radome design based on the measured and computed results. 
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II. BACKGROUND 


A radome is designed to protect the antenna, but not to 
interfere with its operation. Radome performance thus depends 
on the material used in the construction. Initial radome 
design in the 1940's incorporated plywood as the radome 
material. However, this material proved to be unsatisfactory 
due to its high moisture absorption in wet weather and the 
resulting degradation in electronic transmissions. Fabric 
reinforced resins were then developed and exhibited good 
electrical properties. The one drawback to using resins was 
their poor resistance to physical stress. In order to 
increase the strength-to-weight ratio early designers simply 
increased the thickness of the material. However, the 
importance of maintaining thin radome walls with respect to 
the radome wavelength to minimize radome wall reflections was 
not appreciated at the time. This was clearly illustrated by 
the U.S. "pathfinder" program. The poor performance of a 
radar equipped bomber to guide bomber formations to the target 
almost doomed the progreim in its genesis. Problems such as 
this clearly indicated that there was a need for a rational 
and generalized theory of electrical design for future 
radomes. [Ref. 4; pp. 1-3] 

Today there is an even greater need for an understanding 
of a radome and its ability to act as a window for an antenna. 
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The advent of radar guided missiles meant that radomes had to 
meet the physical stresses of high speed flight, while at the 
same time providing optimal electrical performance. Nowhere 
is this more important than in a high speed anti-radiation 
missile, which is designed to seek out and destroy radar 
antennas by homing in on their emitted RF energy. 

A. HISTORY OF ANTI-RADIATION MISSILES 

Historically, the first effective use of an anti-radiation 
missile (ARM) occurred during the Vietnam conflict. The two 
missiles that were widely used during this period were the 
AGM-45 Shrike and the AGM-78 Standard ARM. These two missiles 
had only limited success due to their narrow frequency 
coverage and enemy tactical countermeasures. Efforts to 
improve these systems resulted in the High Speed 
Anti-radiation Missile (HARM). 

1. THE AGM-45 SHRIKE 

The Shrike (Figure 2.1) is a medium range first 
generation ARM that entered service in the mid - 1960s. Its 
development was directly in response to the growing radar 
guided missile threat in Vietnam. On board aircraft ESM/ECM 
receivers provide initial warning of illumination by ground 
defense radars. Once acquisition of the radar is achieved by 
the missile warhead it is then fired. Guidance is 
accomplished by passive radar homing using one of twelve 
different homing head assemblies set to a specific frequency 
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range. It was highly successful in Vietnam, the Middle East, 
and in the U.S. Navy Libya strikes in 1986. 



Figure 2.1. AGM-45 Shrike. 


2. THE AGM-78 STANDAIU) ASM 

The Standard ARM (Figure 2.2) was developed in the 
late 1960s and is an air launched version of the RIM-66 
Standard surface-to-air missile. It has increased range and 
a larger warhead than Shrike, and contains a single homing 
assembly programmed by the aircraft operator. 



Figure 2.2. AGM-78 Standard ARM. 


3. THE AGH-88 HARM 

The HARM (Figure 2.3) is a second generation ARM 
designed to replace both the Shrike and the Standard ARM. The 
main improvements incorporated included a broadband antenna, 
higher speed, programmable software, and increased frequency 
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range. Unlike Shrike, HARM uses a single warhead with a 
broadband capability. This allows the missile to operate 
against a large number of RF threats. Its software controlled 
processing gives it the cd)ility to be reprogrammed in the 
field as opposed to using hardware modifications in the 
missile as enemy weapons develop. The seeker head employs a 
fixed antenna that greatly simplifies the design and cost. 
The seeker head is comprised of the antenna array, ten 
microwave circuit boards and a video processor. This 
configuration provides direction finding capability without 
the need for a moving gimballed assembly. The antenna array 
consists of a plane spiral helix with other elements within 
a radome which itself forms an integral part of the overall RF 
antenna system. This missile is capable of being carried by 
a wide range of aircraft and is used by both the U.S. Navy and 
Air Force. [Ref. 5] 



Hlllllllllllll 
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mm 






Figure 2.3. AGM-88 HARM. 
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B. THE ANTI-RADIATION SEEKER 

The ARM seeker passively detects the power emitted from 
the victim radar. The power received at the antenna is given 
by the radar equation 

- (-M 

P, is the transmitter power, R is the range from the target 
transmitter to the anti-radiation missile seeker receiver, and 
G, and G, are the gains of the transmitter and receiver 
antennas. In practice, the radome introduces an additional 
loss in gain relative to the antenna in free space. [Ref. 3: 
p. 18] 

The seeker design is both challenging and difficult as it 
must cover an extremely wide frequency range. ARM seekers must 
also have the ability to lock on to both the side and back 
lobes of the victim radar. The ARM is guided to the proximity 
of the victim radar by guidance signals generated by the 
antenna, processing circuitry, and computational equipment. 
The antenna is a crucial part of the seeker system because it 
provides the radar signals which are then processed to 
determine angle of arrival. However, the antenna information 
can potentially be distorted by the radome that protects it. 
Only a few radar signal parameters are used by the seeker to 
sort and classify the threat pulse radars by type. [Ref 6; pp. 
18-24] 
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III. RADOME DESIGN 


Broadband, multi-octave frequency coverage is a 
challenging requirement in the design of the modem ARM 
seeker. It allows the seeker the flexibility to scan all the 
frequencies that the target radar will radiate. As a result, 
the radome of the missile must be able to operate over the 
same frequency range of the target emitters. The transmission 
requirements for a radome on a passive radiation sensing 
missile are not as demanding as compared to the requirements 
for an active or semi-active radiation sensing missile for 
equal range performance and moderate power transmitter 
capabilities. Unlike the active or semi-active systems the 
passive seeker on an anti-radiation missile operates on a much 
higher signal-to-noise ratio. For the active and semi-active 
systems, if maximum range is the design goal then the radome 
must pass at least 90% of the incident energy. This is due 
primarily to the low signal-to-noise ratio that these seekers 
must work with. The passive sensor, however, can operate with 
a radome that has up to 50% transmission loss. [Ref. 7; p.l6] 

A. AERODYNAMIC CONSIDERATIONS 

Radome design is generally dictated by aerodynamic 
considerations. Reduction of drag is of prime importance when 
considering speeds more than Mach 1. The conposition of the 
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radome material plays a key role in determining aerodynamic 
heating, dynamic pressure, and rain erosion. The ideal 
electrical shape for a missile is a hemisphere, however a long 
and slender shape is the best aerodynamically for a missile. 
A compromise must therefore be reached in which the 
aerodynamic shape does not degrade the electrical performance 
of the radome. [Ref. 7; pp. 1-4] 

B. RADOME ELECTRICAL CHARACTERISTICS 

RF energy transmitted through the radome is a function of 
shape, frequency, wall construction, incident angle, and 
polarization. Given aerodynamic constraints, the only 
parameter that can change in the design is wall construction. 
The number of layers, the dielectric constant of each layer, 
and the thickness must be taken into account in the design of 
the radome wall. 

Single material radomes are the most common radomes found. 
The thickness of these radomes is adjusted so that it is close 
to one-half the wavelength of the electromagnetic wave 
propagating in the material at the frequency of interest. 
This condition minimizes multiple reflections on the inside of 
the material. These radomes provide good phase 
characteristics and high transparency. However, at low 
frequencies X/2 may be too thick to be practical, whereas at 
high frequencies it may be too thin. 
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incident wave 



reflected wave 

radome 

wall 


impedance 

looking 

z 

in 


transmitted wave 

medium 1 (free space) 

•no 

medium 2 

•n. 

medium 3 (free space) 

' 0 


Figure 3.1 Single boundary layer. 


A crude estimate of the reflection loss (Figure 3.1) due 
to the radome is given by the plane wave reflection- 
coefficient at an air/dielectric interface. The intrinsic 
impedance of media 1 and 3 is r/f, and the intrinsic iit^edance 
of medium 2 is j/j. The impedance looking into the radome is 

denor-vd The reflection coefficient seen by an incident 


wave is 


P = 


^ia~^0 


(3.1) 
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Large reflections are produced when is significantly 
different than ij, . All the energy is transmitted when Z*, is 
equal to ij#. 

In practice the reflection from the second boundary is 
important and has significant effect on the value of Z*,. For 
a given incidence direction, frequency can be chosen so that 
the reflection from the second surface will cancel with the 
reflection from the first surface. This is referred to as a 
tunned radome, but this method is only effective over a narrow 
frequency range. An example of transmission through a tunned 
radome is shown in Figure 3.2. 



Figure 3.2 Transmission through a tunned radome. 
Dielectric constant of material e,= 3.5, reflection 
coefficient 6 = .05, and incidence angle ^(=0. 








Using impedance matching, broadband transmission can be 
achieved for the radome. The outside layer of the radc^ne is 
constructed of a material with a low dielectric constant. Each 
subsequent layer has a higher dielectric constant them the 
previous layer. This creates essentially a tapered 
transmission line as illustrated in Figure 3.3.[Ref. 7: pp. 
16-18] 



Figure 3.3 Stepped inpedance transformer. [From Ref. 7] 

Another technique in designing the wall of the radome is 
to alternate the layers of low and high dielectric constant 
material as shown in Figure 3.4. There are several types of 
laminated walls that enploy this technique of alternate 
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Two - Ply A • Sandwich 

Sandwich 



B - Sandwich 

Multilayer 


# Porous Material 

C - Sandwich 

(_) Dense Material 


Figure 3.4 Lcimlnate radome material. 


layers. The first is the two-ply sandwich which is composed 
of a high dielectric material on the outside with a low 
dielectric material on the inside. However, this simplest 
form of lamination is mainly advantageous at low frequencies. 

One of the most frequently used multi-layer designs is the 
A-sandwich. It is composed of a low dielectric core with a 
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thin skin of denser material on both sides. It provides a 
rigid structure with good strength-to-weight ratio and 
broadband capcUoilities. 

The B-sandwich is a single layer of high dielectric 
material sandwiched between two low dielectric materials. A 
design with the weaker low dielectric material on the outside 
is generally not used. The low tolerance of the exterior 
material to in-flight stress and environmental conditions are 
the main drawbacks despite its being superior to the A 
sandwich electrically. 

The C-sandwich consists of multiple layers of low and high 
dielectric material. This design is inherently structurally 
rigid and is used where broadband capabilities, as in an ARM 
radome, are the prime consideration. Plastic laminates are 
used for the different layers providing the high degree of 
strength required. [Ref. 4: pp. 44-47] 

The radome to be used for the gain and pattern comparison 
is an AGM-88 HARM radome. This radome has a slightly modified 
ogive shape with an A-sandwich wall construction having an 
the inner core made from a soft honeycomb. This con^osition 
of the radome wall offers broadband capabilities to the 
seeker. The nominal thickness and dielectric constants of the 
HARM radome are illustrated in the Figure 3.5. [Ref. 8: pp. 1- 
5] 
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1. DBTEKMININ6 MDLTILAYER RADOMB HIPBDANCE 


£=3.5 1.2 3.5 

r 


exterior 

1 



interior 

5 


Width (mm) 1.37 3 0.7 


Figure 3.5 HARM radome wall construction. 


From Figure 3.5 it can be seen that each layer of the 
radome has a different dielectric constant. The transmission 
line equivalent illustrated in Figure 3.3 can be applied to 
this multi-layer radome. Let rif be the intrinsic impedance of 
the region i as defined in Figure 3.5, and Zf be the impedance 
looking into the region i (at the boundary between the regions 
i and i+1) . Using this notation, the impedance ^ for normal 
incedence is 




+ j'i1_^tan (pi tj) 1 
(Pj ti) 


2 ^4 


(3.2) 
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In equation (3.2) 


= 



(3.3) 


Pi = « v/^o€o€r. 


(3.4) 


and Cj is the thickness of section i. 

The radome interior is region 5 and the exterior 
region 1. The input impedance seen by a wave from the 
exterior region is computed by transforming the "free space" 
load of the interior back to the exterior surface. An input 
impedance is computed at each interface and then used as a 
load to be transferred back to the previous interface. This 
is repeated until ^ is computed for the final outer layer. 
The impedance looking in 2^ (2^ «= in the notation of 
equation (3.2)), is then used to solve for the reflection 
coefficient using equation (3.1). Note that an equivalent 
single layer radome could be determined from the equation 





tlo-^JTletan (P^tJ 
TJ^ + JTlotan (Pe^^e) 


(3.5) 


Since is con^lex there is enough information to solve for 
both tf, and /S,t,. [Ref. 9: pp. 454-467] 


C. OGIVE 6EOOCEETRY 

The most common shape used in the design of missile 
radomes today is the ogive shown in Figure 3.6. The ease of 
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fabrication and acceptcUDle high speed performance make this a 
popular shape. The ogive is a segment of a circle with a 
radius {R) larger than the base radius (b) of the radome. The 
curve generated for the radome is rotated around the missile 
axis (z) forming a body of revolution (BOR) . The ogive shape 
is given by the equation 

r = y/R^ -z^ b - R . (3.2) 


[Ref. 4: pp. 47-48] 



Figure 3.6 Ogive geometry. [After Ref. 1] 


The curve of the HARM radome analyzed in this thesis is 
not a pure ogive, but a composition of four different curves 
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as shown In Figure 3.7. The radius r is a function of the 
coordinate z as follows: 


Segment A; 

= (.7002z) + . 13 

(3.3) 

Segment B: 

Zg = 1.312/ri^'(n^57 

(3.4) 


Segment C: r^=v'( 21 .7964) ^ - {^-11.3988)^ ( 3 , 5 ) 

- 18.1152 


Segment D: (213. 1228 ) ^ - (z-36.1048)2 (3 gj 

- 207.8924 



Figure 3.7 HARM radome segments, where 1 unit is 
approximately 2.05 cm. 
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As shown in Figure 3.8 an ogive with 12 = 17 closely 
matches the actual HARM radome curve. This approximation 
will be used in the computer simulation for convenience. 
[Ref.8: pp. 1-5] 



Figure 3.8 Plot of the ogive curve with R = 17 and the HARM 
ogive curve using equations (3.3 - 3.6). 

1 unit = 2.05 cm. 
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ZV. PATTEKH, GAIN AMD VSNR MEASUREMEMT 


A. PATTERN AND GAIN MEASUREMEMT 

Rectangular horn antennas were used for both transmit 
and receive in the test setup. The radome patterns and gain 
measurements were performed in an anechoic chamber. The 
chcunber is an enclosed room whose walls are lined with 
microwave absorbing material to suppress reflections and 
therefore simulate a free space environment. Absorbers are 
shaped in the form of wedges or pyramidal cones to enhance 
absorption of the incident wave through a gradual transition 
from free space into the lossy absorbing material, thus 
providing very low reflection. 

1. TEST SETUP AMD INSTRUMENTATION DESCRIPTION 

An HP X890A was used as the transmitting antenna, and 
a Microline 56X1 as the receiving antenna. The Microline 
56X1 is a standard gain antenna capable of receiving 
frequencies in X band. The dimensions for both antennas are 
given in Figure 4.1. 

The 56X1 was mounted on a rotating pedestal and the 
transmitter tuned to an X band frequency (8.0 to 12.0 GHz). 
The signal received by the rotating test antenna is detected 
and is used to drive a chart recorder which plots the test 
antenna pattern. The "E-plane" is defined by the orientation 
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Figure 4.1 The 56X1 and HP X890A horn antenna dimensions. 


of the electric field vector in the horn aperture as 
identified in Figure 4.1. The H-plane is orthogonal to the 
E-plane. They are also referred to as vertical and horizontal 
planes, respectively, since these are the orientation of the 
electric field with respect to the chamber floor. A diagram 
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of the anechoic chamber test equipment is illustrated in 
Figure 4.2. 























B. VSNR MBASUSEMENT 


VSWR measurements for the 56X1 and the HARM radome were 
done using a HP8510 network analyzer shown in Figure 4.3. The 
network analyzer consists of a sweeper or synthesizer to 
provide the RF signal source. The test set provides signal 
separation and the first frequency conversion stage. The 
HP8510 consists of two instruments, the IF/Detector and the 
Display/Processor. The system components are controlled by a 
dedicated bus allowing the HP8510 to take full advantage of 
the capabilities of the various instruments. [Ref. 10; pp. (3- 
1)-(3-8)] 



Figure 4.3 The HP 8510 Network Analyzer. [From Ref. 10] 
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A test set provides the points at which the device under 
test is connected, and provides signal separation to measure 
the four S-parameters. First the network analyzer was used to 
measure the reflection magnitude and phase of the Isolated 
antenna. Next the radome was installed in front of the 
antenna and the measurement repeated. A con^arlson of the two 
identifies the effect of the radome on the antenna VSWR. The 
measured VSWR for the 56X1 was obtained after first 
calibrating the network analyzer for a 1 port device. VSWR 
measurement is calculated from the return loss by 

p = 10^, (4.1) 

sm = / (4.2) 

(1-p) 

where D is the measured return loss value in dB divided by 20. 
[Ref. 10; pp. (10-1)-(10-4)J 

The manufacturer VSWR specification for the 56X1 is 1.15. 
The VSWR as measured on the network analyzer is shown in 
Figure 4.4 for X-band frequencies. For the antenna with the 
radome in front, the measured VSWR is shown in Figure 4.5. 
The maximxim VSWR is essentially unchanged with the radome 
present. 
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Figure 4.5 Measured VSWR of the 56X1 antenna and HARM 
radome. 
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V. MBTHOD OP AHALYSIS 


The objective of this work is to validate the ccanputer 
code by conqparing cc»nputed and measured patterns and gain data 
for the antenna and radome. The antenna and radome are 
defined in the spherical coordinate system shown in Figure 5.1. 



Figure 5.1 Spherical coordinate system. 
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The pattern can be determined from the total electric 
field in the direction of the observation point • By 
calculating the electric field due to the currents on the 
antenna, and the electric field due to the currents on the 
radome, the total electric field can be determined by suimning 
the two. Changes in the initial current distribution caused 
by scattering back from the radome toward the antenna are 
neglected. The scattering back toward the antenna primarily 
affects the VSWR looking into the antenna terminals. This is 
small as verified by the VSWR measurement in Chapter IV. [Ref. 
2: p. 5] 

Assuming a lossless antenna, the terms gain and 
directivity will be used interchangeaJsly. Used in this 
context, gain is a measure of the ability to concentrate power 
in a particular direction. The term G, is the maximum value 
of the gain and is defined as the ratio of maximum radiation 
intensity to the average radiation intensity 


G = 

. _ _ 




(5.1) 


U(0,<t>) is the radiation intensity and dD = sdji6d6dkt> is the 
differential solid angle.[Ref. 11; pp. 29-43] 

The radiation intensity for an antenna in a particular 
direction is the power radiated per unit solid angle, and is 
proportional to the magnitude of the electric field squared. 
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B is the sum of Jg* (electric field of the antenna) and B (the 
electric field scattered by the radome). Therefore the 
radiation intensity is given by 

Uie,^) =^|£( 0 , 4>|2 . < 5 . 2 ) 

2»1o 

[Ref. 2: p. 6] 

A. COMPUTER MODELING 

The computer model used to confute the antenna and radome 
pattern and gain was developed by Francis [Ref. l] . It is 
based on a method of moments (MM) solution for rotationally 
symmetric radomes that has been modified to handle thin-shell 
dielectric radomes. The model is an extension of the Mautz 
and Harrington body of r'?volution code [Ref. 3] . The original 
code as developed by Francis had several limitations, but its 
capabilities were extended by Klopp [Ref. 2]. A summary of 
the computer code capabilities and limitations is given in 
Table 5.1 and 5.2. 


TABLE 5.1. COMPUTER CODE CAPABILITIES 


CAPABILITIES 

Circular or rectangular 
aperture 

Arbitrary BOR shape 

Arbitrary amplitude and 
phase distribution 

Step approximation to 
inhomogeneous radomes 

All antenna near field 
terms included 

Step approximation to 
variations in thickness 
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TABLX 5.2. COMPOm CODS CJkPABILITIBS 


SBSTRXCTIOHS/IiZMZTATZQHS 


Radome is "thin" 
electrically 

Radome properties 
represented by a con^lex 
surface impedance 

Phase scanning only 
(no gimballed antenna) 

Restricted to moderate 
radome sizes due to 
computer memory limitations 

Ohmic loss inside of the 
radome not included 



1. METHOD OF MOMBHTS 

MM is a numerical procedure for solving the E-field 


integral equation (EFIE) for the unknovm current In this 
section a brief summary of the MM solution will be presented. 
A more detailed derivation can be found in [Ref. 3]. 

To begin with, the radome curve is discretized as 
shown in Figure 5.2. The surface is treated as perfect 





Figure 5.2 Discretized radome. 
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conductor with a correction term to account for the dielectric 
material properties to be added later. On the surface, the 
current basis functions are: 

sZ^. = B n=0,±l,...±« i=l,2. . . ,Ng-2 (5.3) 

Pit) 

j = $— - — -e'j^ n=0,±l, . . .±*>0 i = l,2. . . ,Ng-l • (5.4) 

Tjt) is a triangle function extending over two segments, and 
Pi(t) is a pulse function extending over a single segment. 
These functions are shown in Figure 5.3. The unit vectors B 
and $ are in the radial and azimuthal directions. The 
variable t is the arclength along the generating curve of the 
BOR, and the variable p is the distance of a point from the z 
axis. N, is the number of surface generating points. Thus, 
the electric current jJ on the surface is represented by 




N.-2 


w.-i 


E " E 

p-i gri 


qn 


(5.5) 


The MM testing procedure (Galerkin's Method) reduces 
the EFIE to a matrix equation for the unknown coefficients for 
each azimuthal mode, I, 

la = • < 5 . 6 ) 
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Figure 5.3 Basis functions; pulse (P) and triangle (T). 
[From Ref. 1] 

Zn is a square matrix called the in^edance matrix for mode n, 
which can be partitioned into blocks corresponding to the 
vector directions 6 and $ 










« tt n 

"11 "n 


(5.7) 


The test function direction is indicated by the first 
superscript, and the expansion function direction by the 
second superscript. is a coliiinn vector called the 


excitation vector, 


Its elements are determined from the 


electric field on the surface of the radome, and will be 
discussed in more detail later. Detailed e:q)ressions for the 
elements of Z are given in [Ref. 1]. 

2. THIN SHELL APPROXIMATION 

To account for the fact that the radome material is a 
dielectric rather than a conductor, a correction term ( ) 

is added to the MM impedance matrix 




(5.8) 


Zmm is the MM impedance matrix for the perfectly conducting 
surface. The elements of Z. are determined from 


Fl]!!„ - /JiT, • (p, g. t,, 


(5.9) 
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or, because of the orthogonality of B cuid ^ 



(5.10) 


The functions H, are the MM testing functions used to obtain 
equation (5.7). For Galerkin's method they are the con^lex 
conjugates of the expansion functions iZ,. Due to the choice of 
triangles in t and pulses in 4>, the correction matrix will be 
diagonal for the Z** block of equation (5.7) and tridiagonal 
for the Z* block. 

The radome surface irrpedance is given by 




1 

j« (e-Cp) t 


(5.11) 


where t in this case is the thickness of the shell. It can be 
rewritten in terms of the dielectric constant (e), thickness 
in wavelengths ratio (nX), and the loss tangent (tan 6) 


Z 


S 


_60_ 

nX[j[€j.-l) +ertan6]' 


(5.12) 


This equation is based on a solid homogeneous radome 
material of constant thickness. Multi-layered radomes can be 
modeled as a solid material by using an equivalent dielectric 
constant. This is the approach used to model the three 
layered HARM radome wall. For radomes with variations in 
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thickness and dielectric constant, the surface can be analyzed 
by a stepped index and thickness. On each segment of the 
radome curve, €,, tan 6, and t can be specified independently. 

3. DETERMINING TEE EXCITATICEI VECTOR 

To compute the excitation vector V of (5.6), the 
electric field of the antenna ^must be known, and it can be 
determined by integrating the known current on the antenna. 
The excitation elements are 

'E^ds (p-t,<|>). (5.13) 

For a X polarized uniformly excited aperture lying in the xy 
plane, with a linear phase applied to scan the beam in the 
direction 6, (0,=O) , 

iJ. = . (5.14) 

The following expressions are the Cartesian components 
of the radiated electric field 

, 5 . 15 , 

. e:iki9'cos^sin9,-R) 
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-11.4,. 



E* = 




~~ j^Jipcos^-p'cos^) (psin^-p'sin^O 


(5.16) 


- ^^/J(PC°s*-pW)^ ,5 

• z(32e^*<‘‘'‘=®“*'*^“*'’*’p'dp'd4)' . 

Primed quantities are associated with the antenna; unprimed 
quantities with the observation point (radome) . The distance 
between the two points is 



(p-pO ^ +z2 +4pp/sin2 



(5.18) 


Also used in (5.15) through (5.17) are 


r- - -1-jkR 

" i?3 


^ _ 3+3lkR-k^R^ 

-F- 


(5.19) 

(5.20) 


where k = 2-n/\. 

The rectangular components of the electric field are 
then converted to spherical coirqponents by using the 
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transformation 


E* = £’jfsin0cos<|> + £ySin0sin4> + £/cos0 

= Ejf cos0cos<^+ £’jfcos0sin4)-B/sin0 (5.21) 

E^ = -Ejf sin<j> + £’/cos<|> . 


Finally, the spherical components from equation (5.13) can be 
used to generate the excitation elements by numerically 
integrating equation (5.13) 


viii= 

• [sin(Vj[-0)£«(®+cos(v^-0)^/]dtdi(|) 


(5.22) 


^ dtd^ . (5.23) 

Jo Jti Pi 


It will be assiamed that the antenna radiated field in the rear 
hemisphere (6 > 9(f) is zero. [Ref. 12; pp. 1-6] 
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VI. DATA AHALTSZS 


A. COMPUTER PROGRAMS 

The mathematical theory presented in Chapter V is encoded 
in FORTRAN and run on a Sun SPARCstation™under UNIX. Two 
programs are used to conpute and plot the pattern and gain 
calculations: RADOME.F and 6AIN.F. 

1. RADOKE.P 

The program RADOME.F was originally developed by 
Francis [Ref. l] and called LDBORMM.F. It determined the 
series expansion coefficients of the current induced on 

the surface of the radome, and also the total electric field 
at an observation point. It was modified by Klopp [Ref. 2] 
to take advantage of symmetry allowing the antenna mainbeam to 
scan in both B and 0. 

The program is interactive, prompting the user for 
various radome, antenna and calculation parameters. The 
output from the program is two ASCII files. Outldbor is a 
listing of the specific radome and calculation parameters. It 
includes the input parameters, radome geometry and inpedance, 
and a table of the pattern magnitude and phase. Curcoefsdat 
is a file containing the current expansion coefficients as 
well as geometry data. MATLAB™ formatted files for plotting 
the 6 and ^ components of the electric field are also generated. 
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The subroutine "TAPER" in RADCA^E.F was modified to 
include a cosine cunplitude distribution to approximate the 
actual horn aperture distribution. A quadratic phase error 
was also added in the H-plane as discussed later. 

2. 6AZN.F 

The output file curcoefsdat from RADOME.F is used by 
GAIN.F along with files of gaussian integration constants. 
According to equation (5.1), the electric field is integrated 
over all space to obtain the gain of the system. An output 
file named outlnt is generated which contains the calculation 
parameters and results. [Ref. 2: pp. 36-45] 

B. CCMPUTED AMD MEASURED RESULTS 

The program RADOME.F was used to compute the antenna 
pattern for the 56X1 antenna without the radome. Figures 6.1 
and 6.3 show the computer generated plots for the E and H- 
plane patterns respectively. The horn aperture dimensions are 
1 .8X in the E-plane (xz) and 2.6X in the H-plane (yz). Plots 
of the antenna with the radome were then obtained for both the 
E and H-planes, and are shown in Figures 6.5 and 6.7. Gain 
was computed for all plots and the results are summarized in 
Tables 6.1 and 6.2. 

Measured pattern plots for the antenna are shown in 
Figures 6.2 and 6.4. E and H-plane patterns with the radome 
are shown in Figures 6.6 and 6.8. In all cases f « 10 GHz. 
From equation (3.2), the radome impedance is Z = 218 - jl33. 
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Figure 6.1 Calculated pattern plot of the 56X1 antenna in 
the E-plane. Vertical scale in dB and the horizontal scale 
in degrees. 



Figure 6.2 Measured pattern of the 56X1 antenna in the E- 
plane. Vertical scale in dB and the horizontal scale in 
degrees. 
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Figure 6.3 Calculated pattern of the 56X1 antenna in the H- 
plane. 
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Figure 6.5 Calculated pattern of the HARM radome and 56X1 
antenna in the E-plane. 



Figure 6.6 Measured pattern of the HARM radome and 56X1 
antenna in the E-plane. 
















Figure 6.7 Calculated pattern of the HARM radome and 56X1 
antenna in the H-plane, 



Figure 6.8 Measured pattern of the HARM radome and 56X1 
antenna in the H-plane. 
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TABLl 6.1 IfXZlILOBB 


MBASUUEMBMTS (S-PLANB) 


QAIM 



56X1 AMTBllNA 

56X1 6 HARM 

RADOME 

E-PLANB 

(COMPUTED) 

17.04 dB 

15.97 dB 

E-plane 
(MEASURED) 

17.00 dB ± 0.3 dB 

16.70 dB ± 0.3 dB 

TABLE 6.2 MAINLOBB 

GAIN MEASUREMENTS (H 

-PLANE) 


56X1 ANTENNA 

56X1 & HARM 1 

RADOME 

H-PLANE 

(COMPUTED) 

17.04 dB 

15.97 dB 

H-plane 

(MEASURED) 

17.00 dB ± 0.3 dB 

16.40 dB ± 0.3 dB 


In the code a rectangular aperture of dimensions a = 1.8X 
and b = 2.6X was used to approximate the horn aperture. The 
amplitude distribution was assumed to be a cosine and a 
quadratic phase error was also added to the H-plane (y 
direction). Thus the aperture field is of the form 

E{x,y) «Jtcos[^) ° (6.1) 
b 

The measured and computed data is generally in good 
agreement. The estimated accuracy of the gain measurement is 
0.3 dB. In general the sidelobe levels and position are 
within a few dB. No background subtraction capability 
currently exists in the system. Some of the sidelobe 
difference can be attributed to chamber reflections. This is 


46 



































evident from the asymmetry in the measured patterns. 


Variations in the impedance had little effect on the 
computed pattern sidelobe levels. The radiation patterns for 
three impedance values are shown in Figure 6.9. 



Figure 6.9 H-plane plot of the HARM radome with three 
different impedances showing little variation in the 
sidelobe gain levels. 
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VII. CQlfCLnSIOiT 


A. SUMMARY OF RESULTS 

The computer generated data for the HARM radome verifies 
that the programs RADOME.F and GAIN.F agree with measured 
results from the anechoic chamber. Typical input values for 
the computer code are listed in J^pendix A. 

The antenna was modelled as a rectangular aperture with 
a cosine amplitude distribution in the H-plane. A quadratic 
phase error was also added in the H-plane to fill in the first 
null and broaden the main beam. As expected, the patterns 
with the radome have higher sidelobes when compared to the 
antenna alone. This is due to the ogive shape which reflects 
energy towards other parts of the radome surface. The 
resulting multiple reflections result in many sidelobes at 
wide angles from the mainlobe of the antenna. In general the 
sidelobe characteristics of the computed and measured patterns 
were similar. 

The gain predictions were also fairly close in agreement 
with measured data. Repeatability and stability was a problem 
with the measured gain because of equipment problems. Also no 
background sxabtraction capability exists and therefore chamber 
reflections were present. Another potential source of error 
is due to boresighting the antenna cuid radome. The alignment 
was simply done "by eye" and is only accurate to a couple of 
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degrees. Also, the data in Table 6.1 assumes that the 
mismatch loss is the saune for the antenna with and without the 


radome. With all of these factors taken into account the 
measured and computed gain data were within 0.5 dB of each 
other, which is about the measurement accuracy of the NPS 
anechoic chaunber. 

B. RECOMMENDATIONS 

Listed below are the following suggestions for further 
research: 

1. Perform additional measurements of other missile radomes 
and compare with the predictions of the computer model; 

2. Vary the antenna type to determine its effect on the 
radiation pattern and code accuracy; 

3. Incorporate a new and more accurate measurement system 
into the Naval Postgraduate School anechoic chamber in order 
to provide more repeatable data; 

4. Extend the study to include other frequency bands; 

5. Investigate other radome shapes and materials. 

This thesis showed that the programs RADOME.F and GAIN.F 
provided accurate estimates for the radiated patterns for the 
HARM radome. As a design tool, it has the potential to reduce 
cost and development time for future radcane designs. 
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APPENDIX A. INPUT DATA FOR RADOKE.F 


A. INPUT DATA FOR 56X1 ANTENNA 

This appendix lists the data file used to connpute the 
56X1 antenna pattern. The progreun RADOME.F is used, but a 
"microscopic" disc is entered as the radome. Thus the 
scattered field is zero and only the antenna pattern results. 


50 









TABLE A.l. INPUT DATA FOR E-PLANE AND H-PLANE 



LETTER CODE FOR FILE 


BOR GEOMETRY 


DISK RADIUS 


NUMBER OF POINTS 


FILENAME IN T (NT) 


FILENAME IN PHI (NPHI) 


FILENAME ANTENNA 


INTEGRATION (NANT) 



PLOTTING INCREMENT 
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B 


mPDT BATA FOR AGM-SS HARM RAOOME 


This appendix lists a typical RADOME.F data file. 

TABLE A.2. INPUT DATA FOR E-PLANE AND H-PLANE 



B-PLANE 

H-PLANE 

LETTER CODE FOR FILE 

c 

d 

BOR GEOMETRY 

1 (OGIVE) 

1 (OGIVE) 

SURFACE CURVATURE 

12.5 

12.5 

ZPRIME 

0.5 

0.5 

BASE RADIUS 

2.65 

2.65 

FILENAME IN T (NT) 

gaus2 

gaus2 

FILENAME IN PHI (NPHI) 

gaus20 

gaus20 

FILENAME ANTENNA 

gausl2 

gausl2 

INTEGRATION (NANT) 



PLOTTING INCREMENT 

1 

1 

HIGHEST MODE 

10 

10 

PHI IN DEGREES 

0 

90 

SCAN ANGLE 

0,0 

0,0 

COMPLEX IMPEDANCE 

(218,-133) 

(218,-133) 

ARAD.BRAD 

.. .. .. ... . . .... 

1.8,2.6 

1.8,2.6 
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APPENDIX B. mPDT DATA FOR GAIN.F 


A. INPUT DATA FOR ANTENNA GAIN 


This appendix lists the input file for the antenna gain 
calculation. 

TABLE B.l. INPUT DATA FOR E-PLANE AND H-PLANE 



E-PLANE 

E-PLANE 

FILENAME IN PHI (NPHI) 

gausl2 

gausl2 

LETTER CODE FOR FILE 

same as used 

same as used 

CURCOEFSDAT 

in RADOME.F 

in RADOME.F 

FILENAME FOR ANTENNA 

gaus20 

gaus20 

INTEGRATION 



THETA START AND STOP ANGLES 

0,180 

0,180 

PHI START AND STOP ANGLES 

0,180 

0,180 
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B 


IHPDT DBTA FOR RADOMB QAZN 


This appendix lists the input file for the antenna and 
radome gain calculation. 

TABLE B.l. INPUT DATA FOR B-PLANE AND B-PLANE 


1 

E-PLANE 

H-PLANE 1 

1 FILENAME IN PHI (NPHI) 

9 aus 20 

gaus20 1 

1 LETTER CODE FOR FILE 

Scune as used 

same as used 

1 CURCOEFSDAT 

in RADOME.F 

in RADOME.F 

FILENAME FOR ANTENNA 

gausl2 

gausl2 

INTEGRATION 



THETA START AND STOP ANGLES 

0,180 

0,180 

PHI START AND STOP ANGLES 

0,180 

0,180 


I 

i 
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APPENDIX C. COOS POS SAOOMS UfPSDAHCX 


!del impedanl.met 
clear 

% confute equivalent radcmie surface Impedcuice 

fl-8. ; 

f2=12. ; 

fstep».l; 

m=0; 

for f»fl:fstep:f2 
m*m+l; 

w»3e8/f/le9; 

%thickness of each region 
dl=.054*.0254; 
d2=.12*.0254; 
d3=.03*.0254; 

% relative dielectric constants 

erl=3.5; 

er2=1.2; 

er3=3.5; 

k=2.*pi/w; 

% assxame lossless radome materials (alpha=alpha2=alpha3=0) . 

Find betas. 

al=0.; 

a2=0.; 

a3=0.; 

bl=k*sqrt(erl); 

b 2 =k*sqrt (er2K' 

b3=k*sqrt(er3); 
gl=al+j*bl; 
g2=a2+j*b2; 
g3=a3+j *b3; 

% characteristic impedance of each section 
z0=377.; 

zl=z0*sqrt(l./erl); 
z2=z0*sqrt{l./er2); 
z3=z0*sqrt(l./er3); 

% impedance looking into each boundary 
zc=z3*(z0+z2*tanh(g3*d3))/(z3+z0*tanh{g3*d3)); 
zb=z2*(zc+z2*tanh(g2*d2))/(z2+zc*tanh(g2*d2)); 
za=zl*{zc+zl*tanh(gl*dl))/(zl+zb*tanh(gl*dl)); 

% reflection coefficient at the input 
%R(m)=((za-zO)/(za+zO)); 

R(m)=za; 

M(m)=abs(R(m)); 

F(m) =f ; 

Re(m)=real (R(m)) ; 

Im(m) =imag (R(m)) ; 

P (m) =atan2 (Iin(m) ,Re (m)) *180./pi; 
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end 

%plot(F,R).grid 
plot(P.M); 
grid 

met a in^edanl.met 

pause 

plot(F,P); 

grid 

meta 

pause 
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